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; ABSTRACT 

(N ■ 

I ' Context. High spectral resolution studies of cool Ap stars reveal conspicuous anomalies of the shape and strength of many absorption 

lines. This is a signature of large atmospheric chemical gradients (chemical stratification) produced by the selective radiative levitation 
I , and gravitational settling of chemical species. 

^ . Aims. Previous observational studies of the chemical stratification in Ap stars were limited to fitting simple parametrized chemical 

' profiles. Here we present a new approach to mapping the vertical chemical structures in stellar atmospheres. 

I ' Methods. We have developed a regularized chemical inversion procedure that uses all information available in high-resolution stellar 

, spectra. The new technique for the first time allowed us to recover chemical profiles without making a priori assumptions about 

, the shape of chemical distributions. We have derived average abundances and applied the vertical inversion procedure to the high- 

^— ■») ' resolution VUT UVES spectra of the weakly magnetic, cool Ap star HD 133792. 

' Results. Our spectroscopic analysis yielded improved estimates of the atmospheric parameters of HD 133792. We show that this star 

I has negligible v^. sin / and the mean magnetic field modulus (6) = 1.1 ± 0.1 kG. We have derived average abundances for 43 ions and 

■"^.^ , obtained vertical distributions of Ca, Si, Mg, Fe, Cr, and Sr. All these elements except Mg show high overabundance in the deep layers 

• and solar or sub-solar composition in the upper atmosphere of HD 133792. In contrast, the Mg abundance increases with height. 

Conclusions. We find that transition from the metal-enhanced to metal-depleted zones typically occurs in a rather narrow range of 
depths in the atmosphere of HD 133792. Based on the derived photospheric abundances, we conclude that HD 133792 belongs to the 
rare group of evolved cool Ap stars, which possesses very large Fe-peak enhancement, but lacks a prominent overabundance of the 
rare-earth elements. 
C/3 . 
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o 



>< 



,1. Introduction ical predictions notwithstanding, the absence of high-resolution, 

hj , high signal-to-noise spectroscopic observations at that time did 

Michaud ( 19701 considered m detail a process of element sep- ^^t permit a direct comparison between observations and dif- 

aration in stellar atmospheres under the mutual action of radia- f^^^^^^ calculations. Later Babel STW^ carried out detailed dif- 

tive acceleration and gravitational settling. He showed that in at- f^^j^n calculations for Ca, Ti, Cr, Mn, and Sr in reff=8500 K 

mospheres with suppressed convection and turbulence this pro- ^j^u^ atmosphere and applied his model to explain an unusual 

cess could lead to chemical anomalies in the form of under- or ^ 3^33 ^ ^.^^^ ^^.^^ ^^^^^ ^^^^^ extremely wide 

overabundances. This radiative diffusion hypothesis was nrst ap- • , , • ^ « » 

,. , . , , , , .c . wings and a very sharp narrow core in the spectra of Ap stars, 

plied to interpret the abundance anomalies inferred for the atmo- d u i' i w • ■ j .u e ■ i • 

^, . ^ ■ , • ,. » /»,.,, Babel s calculations inspired the use of a simple approximation 

^Ph^of magnetic chemically peculiar Ap stars (Michaud et ^^^.^.^^j ^^^^^^^ distribution in the stellar atmosphere in 

al. 119741, which were known to be slow rotators and whose at- f e ^ e ^- t-,. ^ . ^. . j i 

' — ; — ' . n the form of a step function. This parametrization IS now widely 

mospheres are stabilized by the strong magnetic held. , , .• , . ^.vc » j- * i \/ 

^ .7 & fe employed in observational stratmcation studies, for example Mn 

Diffusion processes are not only responsible for the observed ^g^^ ^^^^^ (Sig^j 20OI 1, Si, Ca, Cr, Fe in Ap stars (Wade et 

average atmospheric abundance anomalies but also lead to an in- ^1. 150031 Ryabchikova et al. l2()()2ll2()()5t . 



homogeneous abundance distribution through the stellar atmo- 
sphere depending on the balance between radiative acceleration Despite the encouraging success of the previous attempts 
and surface gravity. Borsenberger et al. ( 198 1 1 calculated Ca and to model chemical gradients using the prescribed parametrized 
Sr abundance profiles for atmospheres with reff>10000 K and shape, real element distributions in the atmospheres of peculiar 
demonstrated the influence of chemical stratification on the pro- stars may differ considerably from the simple step function, at 
file of the resonance Ca 11 3933 A line. These interesting theoret- least for some chemical elements. Some theoretical studies also 

hinted at the possibility of complex vertical abundance distribu- 

Send offprint requests to: O. Kochukhov, e-mail: olegOastro .uu.se tions (e.g. fig. 6 in Borsenberger et al. 198 1 1. Moreover, recent 

* Based on observations collected at the European Southern self-consistent model atmosphere calculations including diffu- 

Observatory, Paranal, Chile (ESO programme No. 68.D-0254) sion (LeBlanc & Monin 2004i> have identified a number of addi- 
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tional effects, like mass loss and weak turbulent mixing, which 
may significantly affect the shape of the vertical chemical pro- 
files but are difficult to predict ab initio. These interesting hydro- 
dynamical effects, treated as free parameters in the current the- 
oretical diffusion modelling, could be potentially constrained by 
comparison with the observed vertical abundance distributions. 

Meanwhile modern spectrographs at large telescopes, such 
UVES at the ESO VLT, have reached remarkable precision and 
spectral coverage, and now allow us to record very high signal- 
to-noise ratio data covering the whole optical region spectra of 
moderately bright Ap stars. This dramatic improvement in the 
quality and quantity of the available observational material sug- 
gests that unprecedented details of the stellar atmospheric prop- 
erties, and, in particular, the vertical abundance distributions, can 
be deduced with the help of detailed modelling of the spectral 
line profiles. These facts stimulated a recent surge of interest in 
the observational analysis of chemical stratification and empha- 
sized the need to improve stratification modelling techniques. 

In this paper we present a new procedure to study verti- 
cal chemical inhomogeneities in Ap-star atmospheres. We have 
chosen HD 133792 (HR 5623, HIP 74181) for our modelling. 
This star was classified as an AOp Sr-Cr-Eu object by Jaschek & 
Jaschek ( 1959l. Mathys Cl990j noted that HD 133792 exhibits 
remarkably sharp lines and no evidence of a strong magnetic 
field. Martinez & Kurtz (■1994'( detected no rapid pulsational 
variability in HD 133792. Despite its brightness and extremely 
low Ve sin /, facilitating detailed spectroscopic studies, the star 
was never analysed with the model atmosphere method. 

The paper is structured as follows. The principles and nu- 
merical details of the vertical inversion procedure are given in 
Sect. |2] Observations of HD 133792 and spectra reduction are 
described in Sect. |3 Fundamental stellar parameters are deter- 
mined in Sect. ^ and the abundance analysis is presented in 
Sect.|5l We derive the vertical distribution of chemical elements 
in Sect. |6] Sect. concludes the paper with the summary and 
discussion. 
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Fig. 1, Depth variation of the iron abundance and opacities for a 
chemically homogeneous (left column) and stratified (right col- 
umn) stellar atmosphere, a) Fe abundance as a function of the 
optical depth in continuum at A 4924 A; b) depth dependence 
of the line centre (filled symbols) and continuous (open symbols) 
opacities. The bottom panels (c) show the maximum error in the 
normalized intensity profile of the Fe ii 4923.93 A line computed 
with different radiative transfer algorithms and different number 
of layers in the model atmosphere. 



2. The vertical inversion procedure 

2.1. Radiative transfer in a cliemicaiiy stratified atmospliere 

The stability of the atmosphere of magnetic stars facilitates ef- 
ficient radiative diffusion and may lead to a buildup of signifi- 
cant chemical abundance gradients in the line-forming regions. 
Theoretical calculations (Babel 1992 LeBlanc & Monin 2004 1 
and interpretation of the observations of cool Ap stars (e.g., 
Ryabchikova et al. 120021 12005 1 point to the presence of steep 
photospheric abundance gradients: variation by up to 2-3 orders 
of magnitude may occur for some elements over a small range 
of optical depths. This extreme vertical chemical non-uniformity 
has to be carefully accounted for in the solution of the radiative 
transfer equation. It is not clear to what extent existing numerical 
radiative transfer schemes are able to cope with a strong depth 
dependence of the line absorption coefficient and to generate suf- 
ficiently accurate theoretical spectra. Since most of the methods 
rely on discretizing the stellar atmosphere into a number of lay- 
ers, performance of the numerical radiative transfer techniques 
is expected to improve with the number of layers. In the context 
of the vertical inverse problem solution attempted here, the best 
radiative transfer algorithm is the one achieving acceptable ac- 
curacy with the least number of vertical zones in the chemically 
stratified model atmosphere. 

We have evaluated the performance of 8 different radiative 
transfer algorithms commonly used in simulating stellar spec- 



tra: the DELO solution with linear (DELOL , Rees et al- IT^g^ . 
quadratic (DELOQ, Socas-Navarro et al. '2000") and Bezier 
spline (DELOB, Piskunov, private communication) interpolation 
formulae for the source function; the Feautrier method (Rees 
et al. ll989> : Hermitian solution (Ruis-Cobo et al. 1999 1; step- 
wise Unno solution (MW, Martin & Wickramasinghe 1979); 
Discontinuous Finite Element method (DEE, Castor et al. 1992) 
and the integral solution of the radiative transfer equation as im- 
plemented in the ATLAS9 model atmosphere code (IS, Kurucz 
1993t. The disk-centre intensity across the Fen 4923.93 A line 
was calculated with each of the algorithms assuming an LTE 
source function and using an identical set of line and continu- 
ous opacities. The latter were calculated using the SYNTH3 code 
(Kochukhov 2007 Ryabchikova et al. l2005l for the 49-layer 
Tes - 7700 K, \ogg - 4.0 chemically stratified model atmo- 
sphere kindly provided by LeBlanc & Monin ( 2004 1. In this the- 
oretical model the Fe distribution changes from log(Fe/A^tot) = 
-3.6 at the bottom of the atmosphere to log(Fe/A^tot) - -9.0 in 
the uppermost layers with a large gradient at logT4924 ~ -1.0. 
For comparison we also evaluated the disk-centre intensity with 
the same model atmosphere and homogeneous Fe abundance 
log(Fe/A^tot) = -4.4. The depth dependence of the iron concen- 
tration and of the line and continuous absorption coefficients are 
illustrated in Fig.[^ and b, respectively. 

Calculation of the Feii 4923.93 A line profile was carried 
out with the original 49-layer depth grid and for the models in- 
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Fig. 2. Normalized disk-centre intensity profile of the Fen 4923.93 A line computed for the stratified iron abundance distribution 
shown in Fig. fusing different radiative transfer algorithms: a) DELOQ, b) DELOB, c) Feautrier, d) MW, e) IS. In each plot the 
lower panel shows the line profiles computed for the model atmosphere with 49 (dotted line), 118 (dashed line), and 332 (solid line) 
layers. The upper panels show the difference with respect to the computation for the model atmosphere discretized into 628 layers. 



terpolated onto a finer grid, containing 118, 332, and 628 lay- 
ers. Interpolation employed the adaptive depth grid refinement 
procedure of the SYNTH3 code, which subdivides layers where 
the line absorption coefficient changes rapidly. The maximum 
error of each of the resulting Fe n line profiles was recorded rel- 
ative to the DELOQ solution for the 628-layer grid. The Fe ii 
line shapes corresponding to the stratified case are presented 
in Fig 12] As expected, little discrepancy is found between the 
spectra computed for the model with the maximum number of 
grid points. This verifies the internal consistency of the radiative 
transfer schemes. However, as illustrated in Fig. [Q;, the algo- 
rithms show rather different convergence properties. In the case 
of the chemically homogeneous atmosphere, most of the meth- 
ods already provide an accuracy of 1 % or better for the original 
49-layer model. At the same time, the stratified case turns out to 
be considerably more challenging: all methods except Kurucz's 
integral solution require at least 100 layers in the model atmo- 
sphere to achieve the necessary accuracy. On the other hand, the 
former method attains an acceptable maximum error of w 0.5 % 
for the original 49-layer model atmosphere. 

The outcome of our numerical tests highlights the stability 
and efficiency of the IS radiative transfer algorithm. In principle, 
each of the studied numerical schemes can be successfully used 
for the forward modelling of the effects of chemical stratification 
if a sufficiently fine vertical discretization is used. The advantage 
of the IS method is to provide the required accuracy at the least 
cost. This is why we consider this algorithm to be the optimum 
choice for the vertical inversion procedure (VIP) which recon- 
structs a non-parametrized chemical distribution using a fixed 
depth grid. Taking these results into account, we have adopted 
the IS algorithm in our vertical inversion code described below. 

2.2. Spectrum s)/nthesis 

Synthetic spectrum calculation for a chemically inhomogeneous 
stellar atmosphere is based on the SYNTHV code written by 
Tsymbal ( 1996 1. The code accepts an LTE plane-parallel model 
atmosphere in the ATLAS9 format. In addition, the input of VIP 
consists of a configuration file that specifies wavelength regions 
of interest and indicates the respective line lists. Atomic line pa- 



rameters are extracted from the VALD database (Kupka et al. 
[T999b . 

The code calculates new number densities for a given 
abundance stratification of any number of chemical elements. 
Partition functions are calculated using the tables given in the 
PFSAHA subroutine of the ATLAS9 code (Kurucz 1993 1. For 
some of the rare-earth elements we have updated these data 
following Cowley & Barisciano ( 1994). For calculations of the 
ionic populations we take into account up to 6 ionization states 
of light elements (atomic number < 28) and up to 4 states for 
heavy elements. Calculation of the continuum opacity is per- 
formed for the blue and red edges of each of the wavelength 
intervals and is based on the ATLAS9 subroutine KAPP (Kurucz 
1993V We include Hi, Hei, Hen, H , different metallic bound- 
free and free-free transitions, Rayleigh scattering for Hi, He i 
and Hj, and electron scattering. 

For each wavelength region the VIP code calculates the line 
opacity due to all atomic lines included in the input line lists. 
Radiative damping. Stark broadening and van der Waals broad- 
ening are taken into account. For all lines, except the He i and 
hydrogen lines, we use the Voigt profile to approximate the line 
opacity and adopt the broadening constants supplied by VALD 
or calculated according to the classical expressions (see Gray 
1992). The Stark broadening of the hydrogen lines is estimated 
by interpolating in the tables by Lemke ( 1997). For the Hei we 
use the Stark broadening profiles of Barnard et al. ( il974t i when 
available. 

The radiative transfer equation is solved with the modified 
ATLAS9 subroutine JOSH. The specific intensity is calculated for 
7 angles between the local outward normal and the line of sight. 
Each intensity profile is convolved with appropriate functions to 
model stellar rotation, radial-tangential macroturbulent broaden- 
ing, and the instrumental profile. The flux spectrum is produced 
by integrating specific intensities over the stellar disk. 

Different spectral intervals are treated independently, which 
permits us to use different instrumental broadening profiles and 
radial velocity corrections. 
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2.3. Inverse problem solution 

We must solve an inverse problem in order to find the distribu- 
tion of a chemical element with depth in the stellar atmosphere. 
We seek the vertical chemical profile that would provide an ad- 
equate description of the spectral line intensities and profiles of 
all spectral features belonging to a given species. This problem is 
similar to abundance Doppler imaging (Kochukhov et al. 200^, 
except that the aim is to obtain the best-fit chemical stratification 
instead of a 2-D picture of the horizontal abundance distribution. 

Similar to Doppler imaging, the vertical inversion belongs 
to the class of ill-posed problems (Tikhonov & Arsenin ri977> . 
A given sample of spectral lines recorded in the optical spec- 
tra of A-type stars is characterized by a very uneven sensitiv- 
ity to the properties of different atmospheric layers. Typically, 
the chemical composition of the line forming regions (optical 
depth -2 < logTsooo ^ 0) is very well constrained, whereas 
very few or no lines sample the uppermost and deepest lay- 
ers. This presents a formidable difficulty for a vertical inver- 
sion algorithm: a multitude of widely different, high-contrast 
solutions can provide a reasonable description of observations. 
Furthermore, inversions are unstable with respect to the initial 
guess, and the properties of solutions, especially the amplitude 
of abundance gradients, will strongly depend on the vertical 
sampling of the stellar model atmosphere adopted in the inver- 
sion. 

A common procedure to alleviate this non-uniqueness and 
instability of the ill-posed problem is to introduce a regulariza- 
tion. This brings in additional a priori information that can be 
used to constrain the problem and decouple solution from the 
vertical sampling used in the model atmosphere. For instance, 
one can look for a solution that maximizes the entropy (Vogt et 
al. ri987.) or impose a requirement that the solution must be a 
smooth function (Goncharskij et al. ll977> . The second criterion, 
also known as the Tikhonov regularization procedure, is appro- 
priate for the vertical inversion problem and is implemented in 
our code. 

Thus, the vertical chemical inversion reduces to the mathe- 
matical problem of finding a depth-dependent chemical distribu- 
tion e(x) that minimizes the function 

i A j 

Here F°^^ and cr(F°^^) are the observed stellar spectrum and its 
uncertainty. F^j"(e) is the theoretical spectrum synthesis for the 
chemical distribution s(x) = log(N/Ntod- We use the logarithm 
of the continuum optical depth at /I = 5000 A as an independent 
vertical variable x = logrsooo. The index / runs over all spectral 
intervals, and A - over the wavelength points in each interval. 
We assign different weights W, to spectral regions according to 
the quality of the respective observational data, or if there is a 
necessity to emphasise or diminish the relative importance of 
fits to the specific spectral features. The second term in Eq. Q 
represents the Tikhonov regularization, and A denotes the regu- 
larization parameter. 

The chemical distribution s(x) is a discrete function that de- 
fines the concentration of a given element in each model at- 
mosphere layer. Its shape is not prescribed in advance but is 
constrained entirely by the available observations and by the 
Tikhonov regularization function. The role of regularization is 
to ensure the stability of the inversion procedure and to provide 
the simplest (smoothest) solution that fits observations. 
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Fig. 3. a) Dependence of the inferred iron vertical distribu- 
tion on the regularization parameter adopted in VIP recon- 
struction. Solid line: reconstruction with the optimal regulariza- 
tion, dashed line: inversion with 20 times higher regularization; 
dashed-dotted line: distribution obtained using 10 times smaller 



regularization. The horizontal dotted line shows homogeneous 
Fe distribution adopted as the initial guess, b) Reconstruction 
of the Fe distribution using optimal regularization and differ- 
ent homogeneous initial guesses (log(Fe/A^tot) between -2.5 and 
-4.0). The solid curve shows the average reconstructed iron 
stratification. The shaded area represents the full range of Fe 
abundance for the inversions started from different initial values. 

We use the Levenberg-Marquardt method (Press et al. ll986t 
as the minimization algorithm. It combines the best features 
of the gradient search with the method of linearizing the fit- 
ting function, which ensures rapid convergence close to the 
minimum. For the VIP code we have adapted the Levenberg- 
Marquardt routine described by Piskunov & Kochukhov ( 200^. 
Trial inversions demonstrate that VIP converges to a stable solu- 
tion in no more than 10-15 iterations. 
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Derivatives of the fitting function with respect to the abun- 
dance in each layer, dT jdej, are required by the Levenberg- 
Marquardt algorithm and are evaluated analytically for the reg- 
ularization function. A numerical approximation is used for the 
first derivatives of the synthetic line profiles: 

dsj Ae 

where we use Ae = 0.01 dex. 

The regularization parameter A is determined empirically. 
We adjust regularization in such a way that i) the solution e{x) is 
smooth but, at the same time, gives a satisfactory description of 
the stellar observations and ii) the solution is independent of the 
initial guess. In practice, these requirements are fulfilled when 
the contribution of regularization to the total discrepancy func- 
tion T is comparable to, but not larger than the first term of 
Eq. Q. 

Figure [^illustrates dependence of the vertical Fe abundance 
inversion in HD 133792 on the choice of regularization strength. 
The vertical distributions recovered with widely different A 
are shown in Fig. |^. The solid curve corresponds to the re- 
sults obtained with the optimum regularization. Regularization 
that is too weak leads to large abundance gradients (dashed- 
dotted curve) and cannot ensure uniqueness of the solution. On 
the other hand, a much stronger regularization (dashed line) 
smoothes the inferred vertical distribution of Fe and does not 
allow us to achieve a low;^^^ in the fit to observations. 

The stability test of the vertical inversion using the optimum 
regularization is presented in Fig. 133. The thick solid line shows 
the average of 16 Fe distributions recovered in the inversions 
started from constant Fe abundance log(Fe/A^tot) in the range 
from -2.5 to -4.0 dex. The dark area represents the full span of 
the individual vertical solutions. Its maximum width is 0.07 dex, 
and the mean width is 0.03 dex. Thus, despite drastically dif- 
ferent initial guesses, regularization has ensured convergence to 
nearly the same Fe stratification. 

In addition to chemical stratification, elements are often dis- 
tributed inhomogeneously over the surface of magnetic chemi- 
cally peculiar stars (e.g., Kochukhov et al. 2004 1. In principle, 
it is possible to compute a synthetic spectrum taking into ac- 
count a non-uniform distribution of elements over the stellar 
surface and chemical stratification at the same time. However, 
such a 3-D abundance mapping procedure requires complete ro- 
tational phase coverage and can only be applied to rapidly ro- 
tating Ap stars. On the other hand, HD 133792 investigated in 
our paper is an extremely slow rotator (see Sect. |3 and lacks 
evidence of spectroscopic variability. Therefore, no information 
can be derived about possible spotted distribution of elements 
over the surface of HD 133792, and this effect is disregarded in 
the present study. 

3. Observations and data reduction 

High-resolution, high signal-to-noise ratio spectra of HD 1 33792 
were obtained with the UVES instrument of the ESO VLT on 26 
February 2002 in the program 68.D-0254. The UVES spectrom- 
eter is described by Dekker et al. (2000 1. The observations were 
carried out using both available dichroic modes. The detailed log 
of the observations is given in Tabled In both the blue arm and 
the red arm the slit width was set to 0.5", for a spectral resolu- 
tion of about 80 000. The slit was oriented along the parallactic 
angle, in order to minimize losses due to atmospheric disper- 
sion. Almost the full wavelength interval from 3030 to 10400 A 



Table 1. Log of UVES observations of HD 133792. 



Date 


UT 


Setting (nm) 


2002-02-26 


07:06:01 


346 


2002-02-26 


07:07:59 


346 


2002-02-26 


06:59:57 


437 


2002-02-26 


07:01:43 


437 


2002-02-26 


07:06:00 


580 


2002-02-26 


07:08:07 


580 


2002-02-26 


06:59:57 


860 


2002-02-26 


07:01:50 


860 



was observed, except for a few gaps, the largest of which are at 
5760-5835 A and 8550-8650 A. In addition, there are several 
small gaps, about 1 nm each, due to the lack of overlap between 
the echelle orders in the 860U setting. 

The UVES data have been reduced with the automatic 
pipeline described in Ballester et al. ( I2000> . For all settings, 
science frames are bias-subtracted and divided by the extracted 
flat-field, except for the 860 nm setting, where the 2-D (pixel- 
to-pixel) flat-fielding is used, in order to better correct for the 
fringing. Because of the high flux of the spectra, we used the 
UVES pipeline average extraction method. 

All spectra were normalized to the continuum with an inter- 
active procedure that employed either a low-degree polynomial 
or a smoothing spline function. 

4. Atmospheric parameters, rotation, and magnetic 
field 

We used the Stromgren photometric indices of HD 133792 to 
obtain an initial estimate of the stellar model atmosphere pa- 
rameters. The observed colours, b - y - 0.026, ci - 0.180, 
m\ - 1.110 (Hauck & Mermilliod 1998 1, were dereddened 
adopting E{B - V) - 0.09. This colour excess follows from 
the reddening maps by Lucke C1978 1 and high-resolution dust 
maps by Schlegel et al. ( 199HJ. Taking into account H/? - 2.866 
(Hauck & MermilHod 1998 1, we have estabhshed T^s = 9334 K 
and \ogg - 3.84 with the calibration by Moon & Dworetsky 
([1985} implemented in the TEMPLOGG code (Rogers 1995 1. 

Model atmospheres for HD 133792 were calculated with the 
ATLAS9 code (Kurucz 1993 1, using the ODE with 3 times the 
solar metallicity and zero microturbulent velocity. With these 
metal-enhanced models we further fine-tuned the stellar parame- 
ters to fit the hydrogen Ha and H/3 lines. This procedure yielded 
the final parameters Tetf = 9400 + 200 K and log g = 3.7 + 0.1. 
The final model atmosphere of HD 133792 has 72 layers and 
covers optical depths from log rgooo - -6.9 to log rgooo = 2.4. 

The projected rotational velocity of HD 133792 was estab- 
lished by synthetic spectrum fitting of the magnetically insensi- 
tive Fe I lines. In particular, for the Fe i 5434.52 A line (mean 
Lande factor z - -0.01) no rotational Doppler broadening ap- 
pears to be necessary after the instrumental smearing corre- 
sponding to the resolution of our UVES observations is ac- 
counted for (see Fig. [4}. The respective upper v'e sin / limit is 
~ l.Okms-i. 

Using the Hipparcos parallax of HD 133792 (jt = 5.87 + 
0.66 mas, PeiTyman et al. 1997 1 and adopting Teff very close 
to the value established above, Kochukhov & Bagnulo (2006} 
determined IgL/Lo = 2.02 + 0.10 and M = 2.80 + O.UMq. 
Comparison with the theoretical stellar interior models (Schaller 
et al. [T9921 suggests that HD 133792 is significantly evolved 
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Fig. 4. Magnetic broadening of lines in the spectrum of HD 133792. In each panel observations (symbols) are compared with the 
theoretical spectra computed for {B}-Q kG (dashed line) and (B) =1.1 kG (solid line). The leftmost panel shows the Fe 1 5434.52 A 
line which has z - -0.10 and, therefore, is almost insensitive to the magnetic field effects. The other panels show Fei, Cri, and Cm 
lines with large eff'ective Lande factors. 



from the ZAMS and is likely to be close to the end of its main 
sequence evolutionary phase. 

One derives R = 3.9 + 0.5 Rq from the aforementioned value 
of the stellar luminosity and T^s - 9400 + 200 K. This radius 
impUes a lower limit of about 200 d for the rotation period if the 
star is viewed equator-on. 

No resolved Zeeman split spectral lines are found in the 
spectra of HD 1 33792, indicating a mean field strength of ^ 2 kG 
(Mathvs .l990 l. Mathys & Lanz (I1992> have also failed to detect 
relative intensification of the Feii 6147 and 6149 A lines and, 
on this basis, have suggested a zero magnetic field for the star. 
However, detailed magnetic spectrum synthesis calculations by 
Takeda ([1991 1 showed this interpretation to be invalid. For mag- 
netic fields weaker than x 2 kG, the diagnostic content of this 
particular Fe ii line pair is highly ambiguous, and, in general, 
there is no correlation between the diff'erence of the equivalent 
widths of the two Fe ii lines and magnetic field intensity. 

With our high-quality spectra we see clear indication of an 
extra broadening, invariably correlating with the magnetic sen- 
sitivity, for many spectral lines. Examples of this are illustrated 
in Fig. |4] for several Fe and Cr lines with large Lande factors. 
We have carried out polarized radiative transfer calculations with 
the SYNTHMAG code (Piskunov 1999 1 to estimate the mean field 
strength. The best fit is achieved for (B) =1.1+0.1 kG, confirm- 
ing our earlier field strength measurement (Ryabchikova et al. 

Recently Kochukhov & Bagnulo ( 20061 detected a marginal 
positive longitudinal field of a; 120 G based on the two low- 
resolution spectropolarimetric observations of HD 133792 with 
FORSl at VLT. The longitudinal field did not change sign nor 
appreciably vary in strength between the two FORS 1 (Bz) field 
measurements separated by 9 months. The large {B} to (B^) ratio 
obtained for HD 133792 hints that currently this star is observed 
at the crossover phase. 

5. Average abundances of HD 133792 

Our line identification is based on the theoretical spectrum cal- 
culated for the whole spectral region 3050-9000 A using the 
line extraction from VALD (Kupka et al. 1999 and references 
therein) and DREAM (Biemont et al. 1999) databases. Atomic 
data on the REE elements compiled in the DREAM database 
were extracted via the VALD interface. Comparison of the 



synthetic and observed spectra allowed us to choose the least 
blended lines for the chemical abundance and stratification anal- 
ysis. 

Spectral lines in HD 133792 are very sharp due to a negligi- 
ble rotation and a rather small magnetic field intensity. Thus, for 
most elements we have measured equivalent widths and carried 
out abundance analysis with the Kurucz's WIDTH9 code, modi- 
fied by one of us (VT) to accept the line lists in the VALD output 
format. Only in the case of either blended lines or lines affected 
by hyperfine and/or isotopic splitting did we employ spectral 
synthesis for abundance determination (Be, Mm, Mnii, Euii). 

The recent experimental oscillator strengths currently in- 
cluded in the VALD database were used for the following ele- 
ments: Tin (Pickering et al. 2001% Mnii (Kling & Griesmann 
2000; Kling et al. 2001), Ndii (Den Hartog et al. 2004 1, Dyii 
(Wickliffe et al. 2000 1. For Ce ii, Ce iii, Nd iii oscillator strengths 
were taken from the DREAM database (Palmeri et al. .20(K?I 
Biemont et al. 2002 Zhang et al. 2002 1, while for Prm the cal- 
culations (Ryabtsev, private communication) based on the ex- 
tended energy levels analysis (Wyart et al. 2006) were used. 
For Euii atomic data are taken from Lawler et al. (2001j|. 
Abundances from the Cr n, Fe ii, Co ii lines were based on the 
oscillator strengths calculated with the orthogonal operator tech- 
nique (Raassen & Uvlings [l998t . Its advantage was discussed by 
Ryabchikova et al. ( 2005t . 

Hyperfine splitting of the Mn i (Blackwell- Whitehead et al. 
2005) and Mn ii (Holt et al. 1 19991 1 lines was taken into account. 
However, its effect was found to be negligible for Mm lines 
(<0.04 dex) and very small for Mm: lines. For the abundance 
analysis we used only Mn ii lines for which experimental oscil- 
lator strengths and hfs constants are available. The Co i lines are 
too weak in the spectrum of HD 133792 to be influenced by hfs. 

The final results of the abundance analysis assuming chem- 
ical homogeneity of the HD 133792 atmosphere are presented 
in Table|2] We have compared them with the atmospheric abun- 
dances in another evolved, slightly cooler Ap star with a very 
weak magnetic field, HD 20441 1 (Ryabchikova et al. l2005l l. and 
in the solar atmosphere (Asplund et al. 2005) . 

The light elements Be, C, N, and O, are underabundant by 
an order of magnitude in HD 133792 when compared with the 
chemical composition of the solar photosphere. The UVES spec- 
tra allowed us to investigate the resonance Be ii lines, which are 
only marginally visible in HD 133792, indicating an order of 
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Table 2. Atmospheric abundances in the Ap star HD 133792 
with the error estimates based on n measured lines. Abundances 
in the solar atmosphere and in the evolved Ap star HD 2044 11 
are given for comparison. 
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magnitude Be deficiency. Thus, on the basis of much more ac- 
curate spectroscopic observations we have confirmed the earlier 
coarse assessment of a Be underabundance in the atmospheres 
of 4 cool magnetic Ap stars (Gerbaldi et al. 1986 1, which con- 
trasts with the very high Be overabundance derived for non- 
magnetic hot HgMn stars (Boesgaard et al. 1982 1. The upper 
limit for the He abundance, obtained from the He i AAAll line, is 
log{He/N(ot) - -2.5, which is 1.5 dex lower than the solar value. 
The large underabundance of light elements is typical of cool Ap 
stars (Gerbaldi et al. il989i Roby & Lambert ll990l . 



Abundances of Mg, Si and Ca derived from the lines of neu- 
trals and first ions differ significantly, which is an indication of 
chemical stratification (see Ryabchikova et al. 2003). 

The Fe-peak elemental abundances are similar in both 
HD 133792 and HD 204411, with more pronounced Fe and, in 
particular, Cr anomalies in HD 133792. This is in agreement 
with the general dependence of the Cr and Fe abundance on 
the effective temperature in Ap stars (Ryabchikova et al. 2004| 
Ryabchikova 2005 1. We did not notice any significant depen- 
dence of the individual line abundances on the wavelength or 
excitation energy for Cr lines. On the other hand, such a de- 
pendence exists for Fe lines: abundances derived from the low 
excitation lines below the Balmer Jump (BJ) are systematically 
smaller than those derived from the lines above the BJ. This is 
another signature of possible stratification of elements. 

Special attention was given to the careful study of the rare- 
earth elements. The evolved Ap star HD 204411 is very Cr-Fe 
overabundant and REE-poor compared to the majority of mag- 
netic Ap stars, where a large REE overabundance is typical. REE 
lines are very weak in the spectrum of HD 133792 and reliable 
abundances can be derived only for the lines of Ce iii, Pr iii, Nd 11, 
Nd III, and Eu 11. We find a REE overabundance of about 1 .5 dex 
relative to the solar composition, with a hint of decreasing over- 
abundance towards the heavier species. 

6. Stratification in 1-10133792 

For the stratification analysis we have chosen chemical elements 
with a large number of unblended spectral lines (Cr and Fe) and 
elements showing the most conspicuous discrepancies between 
the line-by-line abundances derived in the approximation of a 
chemically homogeneous atmosphere (Ca, Mg, Si, Sr). 

The lack of the accurate atomic data precludes a strati- 
fication analysis for some other elements showing discrepant 
abundances. For instance, a significantly different abundance is 
derived from the Con lines with different excitation energies. 
However, hyperfine splitting strongly affects these lines, and the 
hfs parameters are not known for many Co ii transitions. 

Stratification calculations were performed for 6 elements. 
Atomic parameters of the spectral lines employed for each of 
the element are listed in Table|31 The central wavelength, excita- 
tion potential, oscillator strength, and the Stark broadening con- 
stant are provided. For the Cr ii and Sr ii lines marked by aster- 
isks in Table|3 experimental Stark broadening data were taken 
from Rathnore et al. ( 19 84 J and Fleurier et al. ( 1977 ). respec- 
tively. For all other lines the broadening constants were either 
extracted from the VALD database or calculated using the clas- 
sical approximation (Gray 1992 1. 

The final depth-dependent abundances derived for the atmo- 
sphere of HD 133792 using the VIP procedure are plotted as a 
function of logrgooo in Fig.|5] A comparison between observa- 
tions and theoretical spectra is presented in Figs. |51|5] Each of 
these plots shows observations and the best-fit synthetic spec- 
trum computed with the stratified abundance. For comparison 
the best-fit calculation assuming a homogeneous chemical dis- 
tribution is also given. 

Magnesium We have used 4 lines of Mgi, one unblended 
Mg II line and 3 Mg ii doublets to derive the vertical stratification 
of magnesium in HD 133792. This is the only element for which 
abundance increases with height in the stellar atmosphere. The 
concentration of Mg is slightly below solar in the deep layers and 
reaches a 1.7 dex overabundance above logrsooo - -5.0. The 
transition zone extends over ^4 dex in log tsooo. Figure|8]shows 
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Table 3. A list of spectral lines used for the stratification calculations. The columns give the ion identification, central wavelength, 
the excitation potential (in eV) of the lower level, oscillator strength (log gf), the Stark damping constant ("appr" marks lines for 
which the classical approximation was used), and the reference for oscillator strength. 
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3380.707 


2.940 


0.199 


appr 


W 


Cri 


5296.691 


0.983 


-1.400 


-6.120 


MFW 


Srn 


3464.453 


3.040 


0.487 


appr 


W 


Cri 


5297.377 


2.900 


0.167 


-4.307 


MFW 


Srn 


3474.889 


3.040 


-0.460 


appr 


w 


Crn 


5297.606 


10.754 


-0.320 


appr 


RU 


Srn 


4215.519 


0.000 


-0.145 


-5.50* 


w 


Cri 


5298.016 


2.900 


-0.060 


-4.051 


MFW 


Srn 


4305.443 


3.040 


-0.136 


-5.50* 


w 


Cri 


5298.272 


0.983 


-1.150 


-6.117 


MFW 


Sri 


4607.327 


0.000 


0.200 


appr 


LW 


Cri 


5298.494 


2.900 


-0.350 


-3.749 


K93 


Sri 


4811.877 


1.847 


0.190 


appr 


GC 


Crn 


5305.865 


3.827 


-2.160 


-6.599 


RU 















LZ - Lincke & Ziegenbein (T97T}; KP - Kurucz & Peytremann ( 1975 i; AZ ~ Andersen et al. ( 1967 1; SG - Schulz-Gulde 1 1969); G-Ga rz IT9731 
- corrected); BBCB - Berry et al. fT97ll ; AJPP - Artru et al. ( 1981 1; NBS - Wiese et al. 1 1969 1; TB - TO PBASE (Seaton et al.fT992l; BWL - 
Black et al. 1 1972 (; SR - Smith & Raggett (TJSl l; RU - Raassen & Uylings 1 1998 1; K93 - Kurucz jT993l ; M FW - M artin et al. (r988 (; MRW 
-May et al. Il974t ; BPSl - Blackwell et al. ri979l; BPSS - Blackwell et al. ( 1982a (; BPS2 - Blackwell et al. >1982bH; BWLW - O'Brian et al. 
iT99n ; BK - Bard & Kock IT994j ; W - Warner 1.19681 ; LW - Lambert & Warner 1 .1968.) ; GC - Garcia & Campos ( 1988) . 



that an inhomogeneous Mg distribution substantially improves 
the fit to the Mg ii lines, especially the feature at A 448 1 A. The 
mean deviation between observations and spectrum synthesis - a 
parameter that is used to quantify the quality of the fit - is lower 
by a factor of 2 for the model with a stratified Mg distribution. 



Silicon This element shows outstanding stratification signa- 
tures in the spectra of HD 1 33792. We have employed 7 Si n lines 
and one neutral Si transition for the reconstruction of the silicon 
stratification in the atmosphere of HD 133792. We find that be- 
low logTsooo - 0.0 Si is overabundant by 1.5 dex with respect 
to the Sun. On the other hand, the upper atmospheric layers are 
characterized by a 1.7 dex deficiency of Si. The transition zone 
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Fig. 5. Chemical stratification in the atmosphere of HD 133792. 



is extremely narrow and is centred at logT5ooo = -0.2. A dra- 
matic improvement of the fit quality is evident (see Fig. |8} for 
the theoretical spectrum which takes the Si stratification into ac- 
count. Such calculations allow us to reduce the mean deviation 
by as much as a factor of 3 in comparison with the best-fit ho- 
mogeneous model. 

Calcium Ca is another element showing very strong vertical 
inhomogeneity in HD 133792. The Ca distribution infeiTed for 
this star is similar to that of Si, but shows a larger change from 
the lower to upper atmosphere. Ca is overabundant by 1.4 dex 
in the deep layers, whereas the deficiency with respect to so- 
lar composition reaches 2.5 dex in the shallow layers. Thus, the 
abundance of Ca in the atmosphere of HD 133792 changes by 
almost 4 orders of magnitude. We obtain a 3.6 times lower mean 
deviation with the stratified Ca model. These results are based 
on 4 Car and 3 Can lines (see Fig.|9}. The resonance Can K line 
is especially useful in the chemical inversion. 

Chromium Many strong unblended spectral lines of neutral 
and singly ionized Cr are identified in the spectra of HD 133792. 
For the VIP inversion we have employed 9 Cr i and 17 Cr ii lines. 
Very high quality of the UV portion of the UVES spectra of 
HD 133792 permitted analysis of the outstanding Cm 3180.29, 
3421.20 and 3422.73 A lines (see Fig. |7}. We find weak evi- 
dence for the presence of vertical stratification of Cr. This ele- 
ment is strongly enhanced over the whole line-forming region in 
the photosphere of HD 133792. The Cr overabundance changes 
from 2.7 dex in the deeper layers to 2.2 dex in the upper atmo- 
sphere. The stratified Cr distribution leads to a marginal 20% 
of the mean deviation) improvement of the fit to the observed 
line profiles (Fig0. 



Iron The vertical distribution of Fe was obtained from 1 1 Fe i 
and 12 Fe II lines, including the strong line of Fe ii at A 5018.44 A 
(Fig- EJ. Chemical inversion suggests a 1.9 dex enhancement 
in the deep layers and * 1.0 dex overabundance in the upper 
atmosphere, with the smooth transition at logrsooo = -1.0. 
Comparison of the observed and theoretical spectra computed 
for the homogeneous and stratified Fe distribution is illustrated 
in Fig.|6l The case for the Fe stratification in HD 133792 is some- 
what stronger than for Cr: the inferred Fe abundance step reaches 
0.9 dex and the decrease of the mean deviation is more substan- 
tial {a 40%). 

Strontium Our Sr abundance inversion relied on 5 Sr n lines 
and two weak Sri features (Fig.|5). The most important infor- 
mation about the vertical distribution of Sr in HD 133792 comes 
from the shapes of the profiles of strong Srii lines. Accounting 
for the stratification of this element noticeably improves the fit 
to the outer wings of the Srii 4125.52 and 4305.44 A hnes. A 
large overabundance of Sr is derived for the whole atmosphere 
of HD 133792. The Sr enhancement ranges from 5.9 dex deep in 
the atmosphere to 2.5 dex in the upper layers. Remarkably, the 
VIP inversion provides a hint that the Sr vertical distribution is 
more complicated than the step-like profile found for other ele- 
ments. 



7. Summary and discussion 

Thanks to recent advances in the instrumentation at large tele- 
scopes, high-resolution, very high S /N spectra have been ob- 
tained for a large number of Ap stars. This dramatic improve- 
ment of the data quality stimulated a renewal of interest in the 
atmospheric properties of magnetic CP stars. Modern investiga- 
tions of the spectral line intensities and shapes show that many 
peculiar features in the Ap-star spectra originate from an inho- 
mogeneous vertical distribution of chemical elements in stellar 
atmospheres. Thus, magnetic chemically peculiar stars represent 
the only type of non-degenerate stellar objects where direct ob- 
servations and diagnosis of the chemical diffusion signatures be- 
comes possible. Observational study of the chemical stratifica- 
tion in Ap stars is able to provide crucial and unique constraints 
for the theoretical modelling of chemical diffusion in stellar at- 
mospheres and also helps to understand the relation between dif- 
fusion and hydrodynamic mixing processes. 

Until now observational analyses of stratification in peculiar 
stars have been limited to fitting simple parametrized vertical 
chemical profiles to a small number of diagnostic lines (Wade 
et al. 2003 ; Ryabchikova et al. 120021 r2005 1. Most often, chemi- 
cal distributions were approximated with a step function, whose 
shape is described by two abundance values and the position of 
the transition zone. 

In the present study we have developed and applied the first 
assumption-free method of reconstruction of the chemical strat- 
ification in stellar atmospheres. In our technique, the individual 
chemical composition is derived for all atmospheric layers con- 
tributing to the line absorption. Uniqueness and stability of the 
chemical inversion is achieved by applying the Tikhonov regu- 
larization function. This means that our code finds the simplest 
elemental distribution sufficient to fit observations. 

The vertical inversion procedure is successfully applied to 
the weakly magnetic, evolved Ap star HD 133792. We have re- 
constructed vertical distributions of 6 elements. Magnesium is 
found to have an abundance close to the solar one in the lower 
atmosphere and shows an increase of concentration with height. 
All other elements have the opposite vertical distribution: high 
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abundance deep in the photosphere and a lower abundance in the 
upper layers. The transition region is narrow for Si, Ca and Sr, 
and is located close to log T5000 = 0. Fe and Cr are overabundant 
over the whole atmosphere of HD 133792, and the respective 
change of the concentration of these elements is smaller than for 
other species and occurs between logrsooo = -1 and -2. 

Our investigation of the diffusion signatures in HD 133792 
confirms the validity of the step-like, parametrized stratification 
models applied in previous studies of Ap stars. However, we find 
that for some species the transition region can be rather extended 
(Mg), or the vertical distribution is more complex than a simple 
one-step function (Sr). Thus, the applicability of the latter ap- 
proximation is wide, but not universal. In this context, our auto- 
matic regularized inversion approach appears to be more robust 
and should be preferred whenever possible. 



Stratification analysis of HD 133792 is supplemented by a 
detailed chemical abundance study assuming homogeneous ver- 
tical distribution of elements. We have measured the concentra- 
tion of 43 ions of 32 chemical elements. The outstanding char- 
acteristic of the chemical composition of HD 133792 is a very 
large overabundance of the iron-peak elements in its atmosphere, 
combined with a fairly small enhancement of heavy elements, in 
particular REE. Strontium and, possibly, palladium are the only 
heavy elements in HD 133792 showing a very large overabun- 
dance with respect to the solar chemical composition. 

HD 1 33792 is the second star with a detailed abundance anal- 
ysis, including stratification effects, whose position on the H-R 
diagram clearly shows that the star is close to finishing its main 
sequence life. The first such star was HD 204411 (Ryabchikova 
et al. I2005< . Both stars possess weak magnetic fields and have 
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Fig. 9. Same as in Fig.|6]but for Ca (left panel) and Sr (right panel). 



a small overabundance of the REEs compared with many other 
Ap stars in the 8500-9500 K effective temperature region. The 
most remarkable abundance characteristic of both stars is an 
extremely high overabundance of Cr and Fe, in particular in 
HD 133792. 

Cowley & Henry ( |1979> have defined a group of 5 stars char- 
acterized by strong lines of the Fe-peak elements and REE lines 
weaker than in "normal" Ap stars. HD 204411 is a member of 
this group. The H-R diagram position of these objects indicates 
that all of them are close to the end of their main sequence life 
(Kochukhov & Bagnulo 2006 l. 

We were able to extract high-resolution spectra from the 
ELODIE archive' for the three other stars of this group: 73 Dra 
(HD 196502), HD 8441, and HD 216533. Photometric indica- 
tors suggest Teff - 8900-9200 K for these three objects. A pre- 
liminary determination of the Pr and Nd abundance based on the 
strongest lines of the second ions shows that in all three stars 
abundances of these rare-earth elements are comparable to those 
found in HD 133792. In contrast, REE abundances in the young 
star HD 66318 of similar effective temperature, much stronger 
magnetic field (Bagnulo et al. 2003 1 and the same large Cr and 
Fe anomalies are by an order of magnitude higher. For all stars 
these abundance estimates were based on the same lines of Pr iii 
and Nd in. 

Evidently, on the basis of our detailed abundance analysis 
we may extend the definition of the group of stars with weak 
REE lines as a group of evolved Ap stars with small to moderate 



' Ihttp : //atlas . obs-hp . fr/elodie/| 



REE anomalies, and very large Cr and Fe overabundances. All 
stars of this group are lie near maximum in plots of Cr and Fe 
abundance versus Teff (Ryabchikova et al. .2004). Comparison 
between Cr and Fe abundances in evolved Ap stars with weak 
magnetic fields and in the young star HD 66318 with an ex- 
tremely large magnetic field shows that neither evolutionary sta- 
tus nor magnetic field play an important role in creating the Fe- 
peak abundance anomalies. On the other hand, the REE over- 
abundance is clearly less pronounced in the group of evolved 
Ap stars. Further detailed chemical and evolutionary analysis of 
Ap stars in the effective temperature range 8500-10000 K is re- 
quired to understand the role of the magnetic field and/or evolu- 
tion in creating and maintaining REE atmospheric anomalies in 
magnetic peculiar stars. 
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